A polarizable electrostatic model of the N-methylacetamide dimer

Our previously developed polarizable electrostatic model is applied to isolated N-methylacetamide (NMA) and to three hydrogen-bonded configurations of the NMA dimer. Two versions of the model are studied. In the first one (POL1), polarizability along the valence bonds is described by induced bond charge increments, and polarizability perpendicular to the bonds is described by cylindrically isotropic induced atomic dipoles. In the other version (POL2), the induced bond charge increments are replaced by induced atomic dipoles along the bonds. The parameterization is done by fitting to ab initio MP2/6-31++G(d,p) electric potentials. The polarizability parameters are determined by subjecting the NMA molecule to various external electric fields. POL1 turns out to be easier to optimize than POL2. Both models reproduce well the ab initio electric potentials, molecular dipole moments, and molecular polarizability tensors of the monomer and the dimers. Nonpolarizable models are also investigated. The results show that polarization is very important for reproducing the electric potentials of the studied dimers, indicating that this is also the case in hydrogen bonding between peptide groups in proteins. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1933–1943, 2001